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A short-distance heavy quark mass depends on two parameters, the renormalization scale [i con- 
trolling the absorption of ultraviolet fluctuations into the mass, and a scale R controlling the ab- 
sorption of infrared fluctuations. 1/R can be thought of as the radius for perturbative corrections 
that build up the mass beyond its point-like definition in the pole scheme. Treating J? as a variable 
gives a renormalization group equation. We argue that the sign of this anomalous dimension is 
universal: increasing R to add IR modes decreases m(R). The flow improves the stability of con- 
versions between mass schemes, allowing us to avoid large logs and the renormalon. The flow in R 
can be used to study IR renormalons without using bubble chains, and we use it to determine the 
coefficient of the O(Aqcd) renormalon ambiguity of the pole mass with a convergent sum-rule. 



The pole-mass, m po i c , provides a simple definition of 
a mass-parameter in perturbative quantum field theory, 
corresponding to the location of the single particle pole 
in the two-point function. For the electron mass in 
QED TOpoic is used almost exclusively, but for quarks 
in QCD there are two reasons it is impractical. First, 
at high energies, large logs appear which spoil pertur- 
bation theory with m po i . This problem is cured by in- 
troducing the concept of a running-mass m(/i), where 
the renormalizaton group (RG) flow in /i is controlled 
by a mass-anomalous dimension. The second, and more 
serious problem, is that due to confinement there is no 
pole in the quark-propagator in non-perturbative QCD. 
Thus the concept of a quark pole-mass is ambiguous by 
Ara po ie ~ Aqcd ■ This ambiguity appears as a linear sen- 
sitivity to infrared momenta in Feynman diagrams, and 
results in a diverging perturbation series for any observ- 
able expressed in terms of m po ie, with terms ~ 2™n! 
asymptotically for large n. For the heavy quark masses 
(charm, bottom, top) that we study, this behavior is re- 
ferred to as the pole-mass O(Aqcd) renormalon prob- 
lem PJ, where the Borel transform of the series has a 
singularity at u = 1/2. Schemes without this infrared 
problem are known as short-distance masses, and always 
depend on an additional infrared scale R. 

Typically, R is considered as intrinsic to the short- 
distance quark mass definition, mn(fi). Examples are 
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MS : 




R = 




RGI [2] : 


TORGI, 


R = 


mRGi; 


kinetic [3] : 


Wkin, 


R = 


,,kin. 


IS [4] : 


mis, 


R = 


toisCf 


PS [5] : 


mps, 


R = 


// PS 



where Cf = 4/3. Many schemes have R = m, but this 
is not generic. For instance, the lS-mass is defined as 
one- half the mass of the heavy quarkonium 3 Si state in 
perturbation theory, and its R is of order the inverse Bohr 
radius. In the kinetic and the potential subtraction (PS) 
schemes R is set by cutoffs, /J,j ln and /z PS , on integrals 



over a heavy-quark correlator and the heavy-quark static 
potential respectively. Depending on the scales involved 
in a process, schemes with a specific range of /i and R are 
most appropriate to achieve stable perturbative results. 

The goal of this letter is to consider R as a contin- 
uous parameter, and study the RG flow in R of masses 
m(R, n) = rriR(fi). We consider converting between mass 
schemes itia{R, h) and ms(i?',/i) where R <C R' ■ To 
avoid the O(Aqcd) renormalon in fixed-order perturba- 
tion theory a common expansion in a s (/it) must be used, 
which inevitably introduces large logs, \n(R'/R). The 
RGE in R allows mass-scheme conversions to be done 
avoiding both large logs and the renormalon. We show 
this improves the stability of conversions between the 
MS scheme with R = m, and low energy schemes with 
R <C m that are extensively used for high precision deter- 
minations of heavy quark masses || . The solution of this 
RGE is also used to systematically derive a convergent 
series for the normalization of the u = 1/2 singularity in 
the pole-mass Borel transform. 

To start, translate the bare-quark mass in QCD to 
the pole-mass, mbare = Z m m po i e , where UV divergences 
from scales p 2 3> m 2 appear in the mass-renormalization 
constant Z m . The difference between using m po i c and 
any other scheme m(R, fx) corresponds to specifying ad- 
ditional finite subtractions, §m(R,{i). Let 



m-poic = m(R, /i) + Sm(R, (j,) , 

«nfe 

n=l fe=0 



(2) 



5m(R, n) — i^y^ a n k 



4tt 



ln A 



R 



Here a n k are numbers, and a s is in the MS-scheme with 



da s (p) 
din fi 



oo 

= /3M/z)] = -2a s {^)Y,Pn [- 



n=0 



-in 



n+1 



(3) 



We will only consider gauge independent short-distance 
mass schemes for m(R,fi), where 8m eliminates the in- 
frared ambiguity associated to the pole mass. This re- 
quires that <X( n +i)o ~ 2™n! asymptotically for large n. 
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These mass schemes come in two categories. In 
/i-independent schemes (such as RGI, kinetic, IS): 
m(R, fi) = 0. In these schemes an = 0, and 
a n k with k > 1 are determined by the a n = a n o's and 
/3„'s. Thus the a„'s specify the scheme. Masses in the 
other category have a /z-anomalous dimension (like MS), 
and using d/dln/j, m po ie = one finds -t4 m(R,fj,) — 
— Rj fl [a s (fi)]. Here a n \ and a„o are needed to specify the 
scheme and 7^. 7 U does not depend on ln(^/i?), so all 
a„k with k > 2 are determined. For "mass-independent" 
schemes like MS we always have an = 6CV, and a uni- 
versal 7a at leading order (LO). 

Eq. (H) can be used to identify i? for schemes like those 
in Eq. ([j]) . To see that R is related to absorbing IR 
fluctuations into the mass, consider the PS scheme where 

m PS (R)-m pole = -8m PS (R) = ± J ^V(q).(4) 

\g\<R n) 

Here V{q) is the momentum space color singlet static 
potential between infinitely heavy test charges in the 3 
and 3 representations. In m po i — 8m PS the low momen- 
tum part of the potential precisely cancels the infrared 
sensitivity of m po i e , leaving a well-defined short-distance 
mass, m PS . If we increase R from Rq to R\ then 



-8m 



PS 



(Rl) = I dq 




q 2 V(q) 



4tt 2 



dq 



q 2 V(q) 



4tt 2 



(5) 



so additional potential energy is absorbed into m PS , in- 
creasing the range of IR fluctuations included in the PS- 
mass. In other mass-schemes the precise definition of 
R differs, but the interpretation of this scale as an IR- 
cutoff still remains. Another simple example is what we 
call the static-scheme. The static energy, E st&tlc (r) = 
2m po ie + V(r) is free of the 0(Aqcd) renormalon, which 
requires a cancellation of IR sensitivity between the pole- 
mass and the potential energy V(r). Here V(r) is the 
Fourier transform of V(q). A short-distance static mass, 
m stat (i?), can be defined to make this cancellation ex- 
plicit, with Sm stat (R) = -\V(r = l/R). [We use the 
modern definition of the static potential, which becomes 
/^-dependent at O(a^) |J, but does not suffer from in- 
frared divergences. Since O(af) is beyond the order of 
our analysis we refer to the PS and static masses as /i- 
independent. For convenience we also consider the ki- 
netic scheme in the heavy-quark limit, so there are no 
higher powers of R in Eq. (Q).] Since V is attractive, 
increasing R decreases the PS and the static mass. We 
will see that this decrease is universal. It is described by 
an RGE in R where d/dlnRm(R) < at LO in a s (R) 
for all known physical mass schemes. 

The RGE for R. Consider any /x-independent scheme 
where R is a free parameter, such as the PS, static, 
and kinetic schemes. The pole-mass in Eq. (|j) is R- 
independent, so Rd/dRm(R) = —Rd/dR8m(R). To 
avoid having large ln(/x/i?)'s on the RHS we must expand 
in a s {R), Sm(R) = i?£~ =1 a„ [a s {R)/{4n)] This 



yields an RGE for R 



n=0 



a s {R) 



47T 



n+l 



(6) 



For the kinetic scheme the RGE in R was formulated in 
Refs.j| |. To our knowledge the full implications of 
Eq. M) have not yet been studied. We will refer to 7^, 
7?i 72 as the LO, NLO, NNLO, anomalous dimensions 
of the RGE in R. Here 7^ = 01, 7^ = a 2 - 2/3 ai, 



72 



R _ 



03 — 4/?oa2 — 2(i\a\, so they are determined by 
the non- logarithmic terms in Eq. (^). For /i-dependent 
renormalization schemes with a parameter R, the same 
RGE in Eq. (^) is obtained once we set [i = R and de- 
fine m(R) — m(R,R) and Sm(R) — Sm(R,R). These 
schemes have consistent RG flow in the two-dimensional 
R-jj, plane (vertically and along the diagonal). An exam- 
ple of such a scheme is the jet- mass, m) ct (R, fi), defined 
via the position-space jet-function |1C[| . To estimate un- 
certainties in the RGE in R one can set fi = kR in 5m 
and determine 7^(k), then vary about K = 1. 

Since the O(Aqcd) ambiguity in m po i is ^-indepen- 
dent, the derivative of 8m in Eq. (^) ensures that 7_r[q; s ] 
does not contain this renormalon (a key point!). The 
RGE flow in the parameter R takes us from a renormalon 
free mass m(Ro) to the renormalon free mass m{Ri). 
Examples of LO anomalous dimensions are 



( 7( rr b = 4c F , ( 7( f-) =2^, ( 7 £ 

^kinetic = 16CV/3 , ( 7 ^) jct = 2e 7E C F 



4C F , 
(7) 



One can find another suitable scheme from the MS-pole 
mass relation, by taking m(m) — * R in 8m. We call this 
the MSR-scheme. All these 7^ 's are positive. Thus for 
large enough R we have d/dlnRm(R) < 0, and increasing 
R always decreases m(R). This sign appears as a uni- 
versal feature of physical short-distance mass schemes. 
Now, the norm of 7^ does depend on the scheme. For 
a given change AR it determines the amount of IR fluc- 
tuations that are added to m(R). In a different scheme 
an equivalent amount of IR fluctuations can always be 
added to the mass with a different change AR' . To 
see this, consider rescaling R = AR' with a A > 0. 
We demand A ~ 0(1) to avoid large logs. Expanding 
a s (XR') = a s (R') - /3 In Aa 2 (i?')/( 27r ) + ■ ■ ■ gives 



7cf 
, R' 



A7(f, 7* =*bi -2/3070* In A] 



(8) 



7* = A [ 72 fl ~ (4/?o7f + 2/3i7<?) In A + 4/? 2 7t f In 2 A] . 

Thus at LO a scale change in R just modifies the norm 
of 7(f, and we are free to pick A so that 7^ is equal to 
the LO anomalous dimension in some other scheme. The 
condition (7^ /7(f ) ~ 0(1) identifies a class of related 
schemes parameterized by the scale choice for R (our 
A). For a top-quark in the PS, static, and MSR schemes 
we find {7^, 7f , 7^} = {0.348,0. 108,0.231} PS , {0.546, 
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-0.061,0.143} sta \ {0.348,0.213,0.068} MS Vwhere we let 
7^ ee 7 «/(2/3 ) fe+1 and used Refs. g || ||. Here and 
below we use 5 light running flavors. 

To determine the general solution to Eq. (^) first write 



In 



Ri 
Ro 



ai da R 



dtb(t) = G(t )-G(t 1 ), (9) 



where a* = a s (Ri), a R = a s (R), U = -27r/(/3 a J ), and 
< ee —2ir/((3oot R ). For the first few orders 



G(t) = t + 6i In(-t) 



2t 2 



(10) 



where &i = /?i/(2/3 2 ), 6 2 = (/?? - /3 /? 2 )/(4/^), and 6 3 = 
- 2/3 /3i/3 2 + /3 2 /3 3 )/(8/3 6 ). Note that G'(t) = b(t). 
The function G(t) allows us to define 



A QCB =Re G ^=Roe G ^\ 



(11) 



a definition that is valid to any order in perturbation thc- 

(k) 

ory, and corresponds to the familiar definition of Aq CD 

at N fc LL order. Now by making a change of variable we 
obtain the all orders result 



m(Ri) — m(Ro) 



pin Ri 

/ dliJt R lR [a s (R)} (12) 



= Aqcd f °dt 7R (t) 4e" G W 



dt 



where j R (t) = y R [a s (R)]. This is a very convenient for- 
mula for the solution of the RGE in R. 

Lets consider the LL solution. The RGE is 



(13) 



Here j R (t) = — 7o/(2/3oi) and with the leading terms for 
b{t) and G(t), Eq. (ph gives 



a(o) k r t - t 
qcd 7o I .e. 



m(R\) — m(Ro) 



2/3o 

A (0) 
QCD /O 

2/3o 



* — (14) 



;r(o,ii)-r(o,to)] 



Since t\ < to < the integral is convergent. Here and 
below we make use of the incomplete gamma function 



T[c,t] 



dx x 



c-1 



(15) 



and note that differences like the one in Eq. (|14|) have 
no contributions from the cut present in T[c, t] for t < 0. 
To see which perturbative terms the solution in Eq. Jl4| ) 
contains, recall the asymptotic expansion for t — > 00 



rMJ^e- 4 ^- 1 ' 



T(l-c+n) 

n=Q v ' v 7 



(16) 



For c = this plus the LL relation Ag^ D e ( = R yields 



A 



(o) 

QCD 



r[o,t] 



asym 



-2i?^2"n! 

71=0 



P a s (R) 



47T 



n+1 



(17) 



This is a divergent series, but for Eq. (14) we have 
m(Ri) — m(i?o) 



(18) 



7<?#o 



n=0 



2 ^ n^O 



I 2tt 

r /3 ai i»+ 1 
I 2tt 



1 



n\ 1- — > -ln K — 
^ ft Ro 

k=0 



E 

k=n+l 



hi 



ft! Ro 



which is convergent since /3oOt s (Ri) ln(Ri / Rq) / (2n) < 1. 
Eq. (|l8|) displays the problem of large logs in fixed or- 
der perturbation theory for Ri ^> Rq. The RGE in R 
encodes IR physics from the large order behavior of per- 
turbation theory. It causes a rearrangement of the IR 
fluctuations included in the mass in going from Rq to R\ 
without reintroducing a renormalon. 

Using Eq. ([lj) the solution in Eq. ( |l4| ) can be ex- 
tended to include all the terms up to N fc LL. Let 
[t lR {t) b{t) e- G (*) e* (-if*] ee - £~ Sj H)"i, then 



[m{Rx 



t (i? )] NfcLL = A W 



QCD 



E = s 3 {-iy 



e **t>i [r(-S 1 -j,f 1 )-r(-S 1 - J ,i )] 



(19) 



This solution is real. Here AqcdF[c, t] ~ RO(a),~ c ) en- 
codes the suppression of higher order terms, see Eq. ([l6|) . 
Using 7^ ee 7^/(2/3o) fc+1 the first few coefficients are 



So = 7?, S 1 =7f-(6 1 +6 2 )7 R , (20) 

S 2 = 7 2 fl - (61 + 62)7? + [(1+S X )S 2 + (6 2 + 6 3 )/2] 7o fi . 

Connection to m po i c . Eq. ( |l2| ) allows us to study 
the it = 1/2 renormalon in m po i e . This is done by taking 
the limit Rq — > and to = — ln(i?o/AQCD) + ...—> 00, 
where we continue in the upper complex plane around the 
Landau pole. Since \im Ro ^Q m(Ro) = m po \ c , we obtain 



m(Ri 



^polc 



A 



QCD 



d* 7fl (t)|e- fl «. (21) 



Upon expansion in a 8 the RHS reproduces (—8m) in 
Eq. (||). Eq. (|2l]) represents an all order expression for 
the Aqcd renormalon in to po i c • This result can be manip- 
ulated into an equivalent expression for an inverse Borcl 
transform of B(u). We obtain 



m(R) — m l 
B(u) = 2R 

OO 

A/2 z 



poo 

oic = du B(u) e~ u 
Jo 

OO o 



(22) 



" 2 S'"(l-2„).«.- 



^ r(i+&i+*o 



(23) 
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0.5 1.0 1.5 2.0 A, 

FIG. 1: Convergence of the sum-rule for P1/2 for m po i c 



m(R)- m(R ) (GeV) 
-7.0 



_ _ - - LO 




-10.0 



FIG. 2: Top-mass scheme conversion from Ro = 3 GeV to R = 
163 GeV. Shown are fixed order results (LO,NLO,NNLO) and 
RGE results (LL,NLL,NNLL), both in the MSR scheme. 



t(u) = 



^ ff fc (2 M ) fc +Vi(l, 1 + Si + fc, 2 + h-e, l-2u) 



k=0 



(i+h-e)r(k + e) 



Here e G W e -t (_ t )-6i = Y,t=o9l (-*)"', so 5o = 1, 
.9i = &2, .92 = (&| — &3)/2j etc. The normalization Pw 2 



multiplies all terms singular at u = 1/2 in Eq. (22). 
Since 7n(t) is free of the u = 1/2 renormalon, the large 
order behavior of 7^ is dominated by the next pole at 
u = p > 1/2. This implies that asymptotically for large 
fc, 7^ ~ k\ (f3 ) k+1 p~ k . Given that the sum for /3[a a ] 
(and hence Ylidt) converges, Sk ~ 7fe ~ k\(2p)~ k , so 
the sum over fc in P]/ 2 converges. Since Qg(l/2) — 

Y,T=o S k[( 1 + bi-l)T(k + £)]-\ all sums over k are ab- 
solute convergent for u close to 1/2. From Eq. 
the large- n asymptotic behavior for any m(R) — m po ie 



is a 



n+l 



asym 

l n+l 



P1/2 (2/3 ) n+1 Efco 9t r(l+&i-i + n), 
and this series in £ agrees precisely with the behavior ex- 
pected from the Aqcd ambiguity Thus Eq. (23) gives 
a convergent sum-rule for the normalization P1/2. 

P1/2 allows us to test for a u = 1/2 renormalon with- 
out relying on the n/-bubble chain. Any Borel summable 
series of a n 's in Eq. (||) leads to a P 1 / 2 that rapidly goes 
to zero. The largest physical series of a n 's that sums to 
P1/2 = has a tt = 1 renormalon, whereas Pi/% 7^ for 
any u — 1/2 pole. Due to the universality of the u = 1/2 



renormalon of m po i e , its P1/2 is a unique scheme indepen- 
dent number. In Fig. [l] we plot the sum of terms for this 
P1/2 up to fc = (light/blue), k = 1 (medium/green), 
and fc = 2 (dark/red). We show the PS (solid), static 
(dashed), and MSR-schemes (long-dashed), which are 
each generalized to a class of schemes with A <E [1/2,2] 
using Eq. (^). The convergence is clearly visible, and we 
estimate P 1 / 2 = 0.47 ± 0.10. For comparison, the widely 
used light-fermion bubble chain || (large-n / with naive- 
non-Abelianization, nf — > — 3/3q/2), gives an overesti- 
mate, Pi/2 = 0.80. A different series for Pw 2 was derived 

in Refs. |l4j [y|, evaluating (1 - 2u) 1+hl B(u) in an ex- 
pansion about u = 0, at u = 1/2. It gives Pi/ 2 — 0.48, in 
agreement with our result. One can use a^ sym in Eq. (0) 
to define a mass-scheme and study its R dependence [ p^5[ , 
which however suffers from the uncertainty in Pi/2- 

Top-quark mass measurements from jets rely on an un- 
derlying Breit-Wigner, and should be considered as val- 
ues to(Po) m some scheme with R ~ T t 11 . The top 
MS scheme has R ~ 163 GeV » Po so a fixed order con- 
version to MS involves large logs. If we measure the MSR 
mass at Rq and run to R = [m(m)] MS then we directly 
get this MS mass. In Fig. || we compare conversions be- 
tween MSR-schemes with Ro — 3 GeV and P = 163 GeV, 
using a fixed order expansion in ot s (ji) — a s (nR) (dashed 
curves), and the solution of the RGE in Eq. (n9) for 7^ 
obtained with p = kR (solid curves). Varying n gives 
a measure for the residual uncertainty at a given order. 
The plot shows that Eq.([l^) converges rapidly, with flat 
k dependence at NNLL. Also the RGE results display 
better convergence than the standard fixed order expres- 
sions. Comparisons using the PS and static schemes yield 
the same conclusion, with similar convergence. Taking 



the Tevatron mass 172.6 ± 1.4 GeV lf| as m 



MSR 



(Po) 



we obtain m(m) = 163.0 ± 1.3+;^ ± .05 GeV. The first 
error is experimental, the second takes Po = 3^ GeV to 
account for the scheme uncertainty, and the third is the 
uncertainty in our NNLL conversion. This assumes the 
experimental error accounts for hadronic uncertainties. 

The infrared RG analysis performed here can be gener- 
alized to study higher order renormalons and quantities 
other than quark masses. For an infrared sensitive ma- 
trix element of C(Aq CD ) the anomalous dimension will 
have terms R N a™(R), and the corresponding sum-rule 
will provide info on the Borel singularity at u = N/2. 
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